Quantum and classical multiple scattering effects in spin dynamics of cavity polaritons 
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The transport properties of exciton-polaritons are studied with aUowance for their polarization. 
Both classical multiple scattering effects and quantum effects such as weak localization are taken 
into account in the framework of a generalized kinetic equation. The longitudinal-transverse (TE- 
TM) splitting of polariton states which plays role analogous to the spin-orbit splitting in electron 
systems is taken into account. The developed formalism is applied to calculate the particle and 
spin diffusion coefficients of exciton-polaritons, spin relaxation rates and the polarization conversion 
efficiency under the conditions of the optical spin Hall effect. In contrast to the electron systems, 
strong spin splitting does not lead to the antilocalization behavior of the particle diffusion coefficient, 
while quantum corrections to spin diffusion and polarization conversion can be both negative and 
positive depending on the spin splitting value. 

PACS numbers: 73.20.Fz, 72.25. Fc, 71.36.-l-c, 72.25.Rb, 78.35. +c 



I. INTRODUCTION 

Spin dynamics of charge carriers and their complexes 
attracts lately an increasing interest. The issues of spin 
coherence generation, detection and manipulation be- 
came topical during last years li 

The spin properties of cavity polaritons are of special 
interest both from the fundamental point of view and 
due to possible future device applications i^ The quantum 
microcavity is the quantum well embedded between two 
highly reflective Bragg mirrors. In such structures the 
strong coupling between a cavity photon and a quantum 
well exciton takes place, which leads to the formation of 
new quasi-particles: exciton-polaritons. These half-light 
half-matter particles exhibit both photonic and excitonic 
properties. Their polarization (or spin) dynamics is ex- 
tensively studied both experimentally and theoretically, 
see Ref. and references therein. 

The polarization eigenmodes of the quantum microcav- 
ity are so-called TE- and TM-modes where the electric 
or magnetic field vector is oriented perpendicularly to 
the polariton wave vector, respectively. They are split- 
ted by the longitudinal-transverse (also known as TE- 
TM) splitting- which plays a role similar to the spin 
splitting of electron states in quantum wellsi^ It leads 
to D'yakonov-Perel'-like spin relaxation in the collision- 
dominated regime and to spin precession in the collision- 
free regimei^ii One of the brightest manifestations of the 
polariton spin splitting is the polarization conversion or 
the optical spin Hall effect: under the Rayleigh scatter- 
ing of linearly polarized polaritons the scattered particles 
obtain a certain degree of circular polarizationi^ The 
angular distribution of the circular polarization demon- 
strates the second angular harmonics thus reflecting the 
symmetry of the longitudinal-transverse splitting. 

Spin splitting can strongly modulate the interference 
phenomena inherent to the quantum particles. It is well 
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known that the spin-orbit interaction modifies quantum 
corrections to electron diffusion coefficient or conductiv- 
ity and electron spin relaxation timesi^ i^°'^^ Although the 
quantum interference of excitons has a long history,— li^ 
it became topical only recently with the development of 
microcavitiesJ^i^ 

Here we analyse the spin-dependent interference effects 
in quantum microcavities. We focus on weak localization 
effects in exciton-polariton diffusion and spin diffusion, 
effects of interference in spin relaxation and in polar- 
ization conversion. Our results can be summarized as 
follows: 

1. The quantum correction to the polariton diffu- 
sion coefficient is negative despite the value of the 
longitudinal-transverse splitting. It is in sharp con- 
trast with the case of electrons, where the suf- 
ficiently strong spin-orbit interaction changes the 
sign of quantum correction to the diffusion coeffi- 
cient. 

2. The relaxation of the circular polarization degree 
of exciton-polaritons is enhanced by the quantum 
interference effects while the relaxation of linear 
polarization can either speed up or slow down de- 
pending on the value of the longitudinal-transverse 
splitting. 

3. The efficiency of the polarization conversion can ei- 
ther be increased or decreased by the interference 
effects depending on the value of the longitudinal- 
transverse splitting and relation between the scat- 
tering time and the radiative lifetime of polaritons. 

The paper is organized as follows: in Sec.[TT]we present 
the model based on the kinetic equation for exciton- 
polariton spin density matrix. The quantum corrections 
to the collision integral describing the effects of coherent 
scattering are introduced and calculated in the frame- 
work of Green's function technique. Section HIIl is devoted 
to the calculation of quantum corrections to the particle 
and spin diffusion coefficients of exciton-polaritons. The 



interference effects on cxciton-polariton spin relaxation 
times are discussed in Sec. |IVl The multiple scattering 
effects and quantum interference effects in the polariza- 
tion conversion are discussed in Sec. |Vl 



II. THEORY 

Below we present the kinetic theory of the spin dynam- 
ics of exciton-polaritons with allowance for the interfer- 
ence effects. 



A. Model 

We consider exciton-polaritons formed from the heavy- 
hole quantum well excitons. Their spin projection on the 
structure growth axis z can take values ±1 or ±2. The 
latter states are optically inactive and do not participate 
to the light-matter coupling while the former ones con- 
stitute the radiative doublet. It is convenient to describe 
this doublet as a pseudospin 1/2 statCr^ the pseudospin z 
component describes the emission intensity in the circular 
polarization, and the in-plane components correspond to 
the linear polarization: namely, Sx component is propor- 
tional to the intensity measured in the given axes x — y 
while Sy component corresponds to the intensity mea- 
sured in the axes x' — y' rotated by 7r/4 with respect to 
X — y coordinate frame. 

The (pscudo)spin dynamics of exciton-polaritons is 
most conveniently described within the spin density ma- 
trix approach. It can be represented in a form 

Pk = fk + Sk ■ cr, 

where fk is the particle distribution function and Sk is 
the average spin in the given state k. Under the con- 
ditions of Rayleigh scattering experiments the monoen- 
ergetic distribution of the polaritons is excited and the 
processes of energy relaxation can be neglected.— It means 
that the absolute values of polariton wave vectors fco are 
conserved. We consider the situation where the mean 
free path of the particles / is large enough, kgl 3> 1. In 
this regime the dynamics of scattered particles can be de- 
scribed in the framework of the classical kinetic equation, 
and the quantum effects can be incorporated as correc- 
tions to the collision integralJ^ 

In the steady-state regime the kinetic equations writes: 



h 

To 



Q{fk} 



h Sfc X Ofc + Q{sk} 

To 



9k, 
Qk- 



(1) 
(2) 



Here tq is the polariton lifetime, Q{fk} and Q{sk} are 
the collision integrals, and gk, Qk are the components of 
the generation density matrix 7fc = {gk+Ok'^f) describing 
the particle and spin generation rates, respectively. 



The quantity Jlfc in Eq. ([2]) is the pseudo-spin pre- 
cession frequency related to the longitudinal-transverse 
splitting of polariton modes equal to ?i|r2fc|.— It can be 
written as 



flk = r2(fc)[cos2(^fc,sin2(^fc,0]. 



(3) 



where n{k) is some function of the wave vector absolute 
value k, Lpk is the angular coordinate of fc. It depends 
strongly on the microcavity parameters. Note, that the 
components of Jlfc are described by the second angular 
harmonics of the wave vector angle, because the pseu- 
dospin flip is accompanied with the change of the polari- 
ton spin by two. It is strongly different from the case 
of two-dimensional electrons, where the spin splitting is 
described by the first and third angular harmonics. 

We assume that the scattering of exciton-polaritons is 
caused by a short-range disorder, i.e. the scattering cross- 
section is angular independent. This condition can be vi- 
olated in real structures'^ but our goal is to consider the 
simplest case which allows an analytical solution. Fur- 
thermore, we assume that the polaritons are described 
by a parabolic dispersion, Ek = h?k'^ /2m, with an effec- 
tive mass 771, and introduce the density of states per spin 
V = m/2TTfi^. 

In typical microcavities under the conditions of 
Rayleigh scattering the kinetic energy of polaritons, Ek , 
is of the order of several meV, while the longitudinal- 
transverse splitting is of the order of tenths of meV. 
Therefore the effect of the polarization on the orbital 
dynamics of polaritons can be neglected. At the same 
time, the spin precession frequency, inverse lifetime and 
the scattering rate can be comparable. 

The collision integrals entering Eqs. ((T|), ^ can be 
written as a sum of the classical contribution 

Qcdfk} = Qj2(h-fk')KEk-Ek'), (4) 

k' 

Qci{sk} = Q^{sk- Sk')S{Ek- Ek'), 

k' 

and the quantum corrections. Here Q is the elastic scat- 
tering constant, QT) = t^ with ti being the momentum 
scattering time. 

The quantum corrections to the collision integrals can 
be most conveniently found in Green's function tech- 
nique. Various contributions to the scattering cross- 
sections are exemplified in Fig.[TJ Solid lines are retarded 
and advanced exciton-polariton Green's functions which, 
with allowance for polariton spin and the longitudinal- 
transverse splitting have a form of 2 x 2 matrices and 
readi^ 

5^'^(fc,tj) = [hw-Ek-n{cT-nk)/2±\h/2T]-^, (5) 

where 

Note that the radiative lifetime of exciton-polaritons 
plays role of the phase relaxation time in the theory of 
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FIG. 1: Examples of irreducible diagrams which contribute 
to the collision integral, (a) Single scattering, (b) coherent 
backscattering, (c) and (d) coherent scattering by an arbitrary 
angle. 



electron weak localization^^ The dashed line in Fig. [T] is 
the scattering amplitude which reads hQ/2iT ~ h? /mTi. 
The diagram Fig. [Ija) shows the single scattering pro- 
cess. The diagram in Fig. [Ijb) describes the interference 
of polariton which passes the same configuration of scat- 
terers in the clockwise and counterclockwise directions 
and propagates exactly backwards. Diagrams Fig [Ijc) 
and (d) describe the same interference but accompanied 
by the scattering by an arbitrary angle. The corrections 
to the collision integral can be expressed in the terms of 
the Cooperon operator Cg^{q) which is the sum of di- 
agrams depicted in Fig. [IJb) with any number of lines 
N ^ 3.(^) The Cooperon satisfies the equation 

/3'7' 



with 



hQ 



Vt,\q)='^Y.^%{k',E,)gl{q^k\E,). (7) 



27r 



Here it is assumed that the polaritons have the same en- 
ergy Eq determined by the excitation. As a result the 
distribution functions can be reduced to the angular de- 
pendent parts only 

fk = f{Vk)6{Ek - Eo), 
Sk = s{ipk)S{Ek ~ Eq). 

The shape of the coherent backscattering cone ob- 
tained by the summation of the diagrams in Fig. [H^b) 
is given in the limit of Qt = by the following expres- 

1 2 

sioni^ 



/i(fc)=/°(fc) 



To V(n/r)2 + (fc + fc^)2;2 _ I 



(8) 



Here /°(fe) describes the classical intensity distribution. 
In the case of exact backscattering fc = — fco and tq 3> ti 



^ The diagram with two intersecting impurity Unos is dominated 
by the pairs of impurities separated by ~ fcg^ and can not be 
treated within our approximation kol 2> 1. 



the coherent backscattering intensity is I^{k) — /"(fc). If 
To ^ Ti the coherent backscattering is negligible /^(fc) <C 
/"(fc). The backscattering cone angular width is small, 
of the order of (fco/)^^ ^ 1 and the details of its shape 
are irrelevant for our consideration. Thus, the coherent 
backscattering processes simply couple states with k and 
-fc. 

The quantum effects are most pronounced in the mul- 
tiple scattering regime where tq '^ t. This assumption 
is used hereafter in the description of the interference 
phenomena. Furthermore, in the limit ln(To/Ti) ^ 1 the 
sum of the diagrams Fig.[ljc) and (d) weakly depends on 
the scattering angle. Even if this condition is violated the 
classical effects of finite correlation length of the potential 
leading to the angular dependent scattering cross-section 
may be more important. Therefore their angular depen- 
dence is neglected, and the quantum corrections to the 
collision integral can be written asiiiii 

Q,nt{fk} = -WoJ^ [/(^ - tt) - fiv'MEk - /?o), (9) 

Qyndsk} = -^ j%\<^ -^)- siv'MEk - Eo). 

Here the quantities Wq, W describe spin-dependent re- 
turn probabilitiesiS. 
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where ct* for i = x,y, z are the Pauli matrices. 

We consider an isotropic spin splitting, i.e. Q^. is inde- 
pendent of the angles of the wave vector fc and equals to 
ri(fc). see Eq. Q. Thus, for the cylindrical symmetry of 
the problem under study, the only non-zero components 
of W arc 



W^^W^ 



Wy 



Wn = W, 



Thus, the decription of the spin dynamics of the exciton- 
polaritons is reduced to the solution of the kinetic equa- 
tions with the collision terms in the form of Eqs. ([4]), ([9]). 
Latter depend on the spin-dependent return probabilities 
Wo, W|| and W± which can be found straightforwardly 
from the Cooperon operator C"^{q) Eq. ©. 

B. Cooperon 

In order to find the Cooperon we follow the procedure 
outlined in Refs. l2Clll2ll and make use of the fact that the 
operator V can be presented as 



2tt 



V ^ — — — [1 - itL ■ ftk + iTVk ■ q] 
Ti J It: 




(11) 



where 



C. Spin-dependent return probabilities 



■'Oi'^ .fib 



'■a/3 — CTb-i 



is an operator of the difference of spins of two inter- 
fering particles and vy, is the velocity operator. This 
result can be contrasted with the situation realized for 
two-dimensional electrons where the spin splitting is odd 
function of fe and, thus, the total spin of interfering par- 
ticles enters into the definition of V . 

Nevertheless, in our treatment we use the representa- 
tion of the total spin of the interfering particles: 07 — > 
SrUs^ where 5 = 0, 1 is the absolute value of the total spin 
S, and ms is its projection onto the z axis (|tos| < S). 
The pair of particles with S" = is in the singlet state 
while S = \ corresponds to the triplet one. 

Since V and, hence, the Cooperon are determined 
by the operator L, two independent contributions to 
Cooperon can be separated. Namely, the matrix V can 
be block-diagonalized, one of the blocks corresponds to 
the pair in the triplet state with rris = and another is 
a 3 X 3 matrix corresponding to two triplet states with 
\ms\ = 1 and a to singlet. The part with (S', nis) = (1, 0) 
reads 



(12) 



v/(ri/r)2 + (g02' 
and the corresponding Cooperon is given by 



Co 



p3 



\-P 



(13) 



The operator V in the basis of three other states, 
(S", ms) = (1, 1); (0, 0); (1, -1), has the following form 



P-S'o ie-^i'^-'i? e-^iv,^^ 
Vx = I ie^^'f'^R P - 25*0 -ie-^^'^-'R 
e'^'^fi'^Si -ie^^'f-R P - So 



(14) 



where ipq is the angular coordinate of q, 
00 
S,n = dx exp {- x) J,n (xql) sin^ i——], (15) 



R = —jn / dxexp (— x) J2 (xgZ) sin (xfir), 


and Jm{y) are the Bessel functions. (^) The Cooperon 
corresponding to these three states is given by 



where / is the 3x3 unit matrix. 



(16) 



^ The integrals JTSj can be calculated analytically but the expres- 
sions are very cumbersome, therefore we leave them in integral 
forms. 



Using Eqs. (jlOp and rewriting the expressions for Wq 
and W^||,± in the same basis as Co, Ci, we have 



Wn 



Wn 



Wi 



— 5:{Tr[£iCi]+Co}, 



koT 



koT 



(17) 



^{Tr[Ci]-Co}, 



— 5]{Tr[£2Ci]+Co}, 



knT 



where £1 is the matrix with the diagonal (1, —1, 1) and 
other elements being zero, and £2 is the matrix with 
unit anti-diagonal elements and other elements being 
zero. Note that Wq can be recast in a conventional 
form as a difference of triplet and singlet contributions 
to C^^iSO^ 

The asymptotic values for the components of W tensor 
can be obtained analytically in the case ln(ro/Ti) 3> 1, 
riri < 1. They read: 



Wo = 

W\ = 



1 



TsA 



To 



r. 



21n^^+ln--ln 

ZlTTkoi \ T T T 

^^f21n^+ln^-ln^ 



(18) 



1 



T, 



'^ll 



To 



2..MV^^^V + 1% 



Here the lifetimes are introduced for spin components 
parallel and perpendicular to the growth axis z: 



1 



1 



1 



02 



1 



— = 0^ + -, = +_. (19) 

Js|| To Js_L ^ To 

In this limit, the spin-dependent return probabilities 
contain logarithmic factors, which is a specific feature of 
two-dimensional systems where the return probability is 
proportional to the logarithm of the ratio of longest and 
shortest allowed travel times. The lower boundary for the 
travel time is obviously ^ t while the upper boundary 
depends on the particle lifetime tq and spin lifetimes T^n, 
Ts_L. The expression for the particle return probability 
Wq differs from the corresponding transport coefficient 
derived for electrons in Ref. |T^. It reflects the specifics 
of the even in fc spin splitting in polariton systems. 

In the regime of Or ^ 1 the leading logarithmic con- 
tributions to W have the following form: 

Wo = W^K -^ (0.O6 + In ^) , (20) 

W\\ « ^ \ , fo.06 - In — 

Clearly, the main (logarithmically large) correction is 
determined by the (S*, m^) = (1,0) contribution to the 
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FIG. 2: (Color online) Components of the tensor W plotted as the functions of the product Q,t, where Q. — Q{ko). The radiative 
lifetime to = lOOri. 



Cooperon. Other states give small additional term. It is 
worth noting that at fir — > oo quantitites Wq , W± have 
the same sign as at f2 — *■ 0, while W\\ changes its sign as 
compared with the case fir ^ 1. 

The dependence of Wq, W\\ and W± on spin sphtting 
value in a wide range of fir variation is presented in Fig. [2] 
at a fixed ratio tq/ti = 100. Panel (a) of Fig. [2] shows 
the non-zero components of W in the case of exciton- 
polaritons (i.e. where the spin splitting contains second 
angular harmonics). The case of electrons where the spin 
splitting contains first harmonics is presented for com- 
parison in Fig.m^b) . Although the overall behavior of W\\ , 
W±, Wq is similar in both cases, there is a strong qual- 
itative difference. Namely, for electrons Wq changes its 
sign with an increase of Sir which is a direct consequence 
of the antilocalization phenomena: the sufficiently strong 
spin-orbit interaction leads to the phase n acquired by an 
electron on the closed trajectories, and the probability 
for electron to avoid its initial point incrcascsi^ The real 
spin of polaritons is 1 therefore the aquircd phase is 27r, 
and the quantity Wq docs not change its sign as a func- 
tion of fir. Mathematically, it is a direct consequence of 
the fc-even spin splitting: the interference is governed by 
the difference spin L of the particles contrary to the case 
of odd spin splitting where the interference is controlled 
by the total spin of the particles S. Therefore, in the case 
of electrons the Cooperon can be separated in the singlet 
and triplet parts with respect to the total spin S which 
enter to Wq with different signs i^ On the contrary, in the 
case of polaritons the part corresponding to the triplet 
with zero projection of the total spin (5, Wg) = (1,0) 
is separated from the Cooperon, therefore Wq keeps its 
sign, while Wh demonstrates the antilocalization behav- 
ior. Therefore, the antilocalization of polaritons does not 
occur. 

Figure [3] shows the dependence of Wq on the ratio 
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FIG. 3: (Color online) The dependence of Wo on the ratio 
T"i/ro calculated for the different values of Q,t. The inset 
shows the dependence of all components of W tensor at fixed 
VLt ^ 1. 



ti/tq calculated for exciton-polaritons. The curves cor- 
respond to the different values of fir. It is seen that 
Wo decreases monotonously with an increase of the ra- 
diative rate. An inset to the figure shows Wq, W\\ and 
WIl calculated as functions of ti/tq at fixed VLt = 1. In 
agreement with Eqs. (fT7)) and (fT5)) only W\\ behaves non- 
monotonously with the decrease of the polariton lifetime. 
Its dependence on ti/tq is determined by the competi- 
tion of two contributions to Cooperon: f2-independent 
Co and fl-dependent Ci. They enter into Eq. (|17p with 
the opposite signs and lead to non-monotonic behavior. 
Cooperons Cq and Ci contribute to Wq and W± with the 
same signs and therefore these components of the tensor 
W are monotonous functions of ti/tq. 



III. QUANTUM CORRECTIONS TO THE 

PARTICLE AND SPIN DIFFUSION 

COEFFICIENTS 

The diffusion coefficient D„ for the quantity n{r, t) 
where n can be either the particle density or one of the 
spin density components s^ , Sy and s^ can be defined 
from the Pick's law 



i?„Vn + j"(r,t)=0, 



(21) 



where j" is the flux density of the corresponding quan- 
tity, r is the coordinate, t is time. 

Let us first consider the particle diffusion coefficient 
Dq. We assume that the distribution function of polari- 
tons can be recast as /(r, fc, t) ~ n{r, t)f{Lpk)5{Ek — ftw) 
and that n{r,t) weakly depends on the coordinate and 
on time. The kinetic equation with the collision integrals 
dm, ® reads 

(/71 71 C/Tl /? A* T) V — 

^^f{vu) + -f{^,) + ^ • -/(^.) + - [/(^.) - f 



nWo J{ipk - 71") - / 



0. (22) 



Here / ~ (27r)~^ /„ f[Lp)dLp is the average number of 
particles. Obviously /(<p) contains only zcroth and first 
harmonics in Lp^. The particle current thus reads 



27r 



r(r,i) = n(r,t)/^^/(^.), 

ZTT TO 



(23) 



and Eq. (fSS]) together with Eq. ((2T|) yields for the parti- 
cle diffusion coefficient 



i?0 = ^c/(l-W^OTl), 



(24) 



where Dd = fi.^fc^Ti/2m^ is the classical diffusion coeffi- 
cient. 

Analogously one can derive the diffusion coefficients 
for the spin components Sz (||) and Sx,y (-L): 



Z?||,i = D,,i{l - VK||,xTi) 



(25) 



From Eqs. ([M]) . ([^5]) one can see that the quantum 
corrections to the particle and spin diffusion coefficients 
are determined by the respective components of the ten- 
sor W. It follows from the previous section that Wq is 
positive for all values of Q.T and ti/tq. Therefore the 
backscattering is enhanced and the quantum interference 
leads to the decrease of the polariton diffusion coefficient, 
cf. (|24p . The same holds for the in-plane spin diffusion 
coefficient, Dj_, which is determined by Wx > 0. The 
quantum correction to the longitudinal spin diffusion co- 
efficient, D||, is determined by W\\ which can be either 
positive or negative depending on Q.T and ti/t. In par- 
ticular, W\\ is negative for large Sir and tq/ti, and quan- 
tum effects will lead to an increase of the longitudinal 



spin diffusion coefficient. For small values of ilr longi- 
tudinal diffusion coefficient is decreased by the quantum 
interference. 

Quantum corrections to diffusion coefficients can be 
extracted from the dependence of these coefficients on ex- 
ternal perturbations which introduce an extra dephasing 
of the particlesii^ These effects can be incorporated into 
the effective lifetime tq. For instance, the temperature 
variation modifies the rates of inelastic processes. Since 
exciton is a neutral particle, the magnetic flux through 
the trajectory is proportional to the electron-hole separa- 
tion. Hence the magnetic field affects the interference at 
Ib ~ flB where Ib is the magnetic length and ae is Bohr 
radiusi^ However, such a field is not classically weak, 
therefore strong diamagnetic effects are dominant, there- 
fore the weak localization corrections to magnetodiffusion 
can be hardy separated. 

In microcavities tq can be efficiently varied with inci- 
dence angle of light proportionally to the photonic frac- 
tion in polariton. It is seen for Fig.[3]that the presence of 
the spin-splitting does not lead to the non-monotonous 
dependence of Wq (and, therefore, of the quantum cor- 
rection to the diffusion constant) on the lifetime. On the 
contrary, the antilocalization behaviour can be observed 



inD 



IV. 



QUANTUM CORRECTIONS TO THE SPIN 
RELAXATION RATES 



The spin dynamics is known to be non-exponential 
with allowance for the quantum interference effects: weak 
localization leads to the appearance of the long-living 
tails in spin polarizationJii^ It is convenient to deter- 
mine the tensor of spin relaxation rates F from the bal- 
ance equation for the total spin in the system at the 
steady-state spin pumping 



(^.- + f) 



(26) 



Here s = (27r)^^ J„ s{(f)Aif is the angular average of the 
spin distribution, and the spin generation is assumed to 
be isotropic 

dk =g5{Ek --Bo). 

We represent the distribution functions separating 
terms of zeroth and first order in quantum corrections: 

/M = rM + /'(¥'), s{^) = s'>{^) + s\^). (27) 

Here the upper index refers to the order in Wq, W . 

In order to solve the kinetic equations (H]) , ([2]) we note 
that the solution of the following equation 



s(^) 



s{v) X n{^) = F(</p), 



(28) 



where ^{ip) = rife writes 



OPTICAL SPIN HALL EFFECT 



[F(^) + rn(^) X F(^) + T'n{^){n{^) ■ F{^))] . (29) 

Neglecting the weak-locahzation effects (i.e. putting 
Wo,Vl^||,Wj_ = 0) one arrives to Eq. (^5]) for the spin 

distribution function s^ with 



F= — 

n 



(30) 



It is enough to consider the cases where g is directed 
along the growth axis of the sample, g || z, or lies in the 
plane of the structure, because tensor T for cylindrical 
symmetry has only two independent components, 



r. 



r. 



describing the longitudinal and transverse relaxation 
rates, respectively. 

If g II z the self-consistency equation for s gives s^ = 
gzTs\\, and hence 



s^{ip) = [g + Tniip)xg]T, 



s\\- 



(31) 



Therefore, the quantum correction for z-component of 
spin reads 



g^T^JnrfW^, 



(32) 



and according to Eq. (j26|) the longitudinal spin relaxation 
rate is given by 



II 



17V (1 + tWi 



(33) 



Since W± is positive the z-component spin relaxation is 
enhanced by the quantum interference effects. This equa- 
tion shows the correction to ri in the longitudinal spin 
relaxation rate has an inverse sign as compared to the 
correction to the spin diffusion cocfhcicnt Eqs. ([25]) . It 
is due to the fact that the spin relaxation rate is gov- 
erned by the relaxation of the second harmonic of the 
spin distribution functioni^ whereas the spin diffusion is 
determined by the relaxation of the first harmonic. 
The calculation of the transverse relaxation rate yields 



n'^T/2 



1 + (17r)2/2 



TWn -TW±{nT)^/2 



1 + (r!r)2/2 



(34) 



This equation shows that at small Qt the quantum cor- 
rections to T± are determined by Wh > and the spin 
relaxation is enhanced by an interference. In contrast, at 
large firi the quantity W\i becomes negative and slows 
dows spin relaxation. Besides, the second contribution 
proportional to W±(flTi)'^ > becomes even more im- 
portant and suppresses spin relaxation as well. 



Under the conditions of the optical spin Hall effect the 
TE- or TM-eigenstate of the quantum microcavity with 
the wave vector fcg is excited and the circular polarization 
in the scattered state k is observedi^ In this case the 
generation rate can be represented as 



9k — gSk,k(n 9k — gSk,ko, 



(35) 



where the Kronecker (5-symbol used here is defined as 



c 27r,^ 

OkM = -:^S[Ek 



Ek„)S{ipk - (fiko)- 



We first consider classical multiple scattering effects in 
the optical spin Hall effect regime and demonstrate that 
the spin relaxation decreases the polarization degree as 
compared with the single scattering regime considered in 
Ref. |5|- Further, we calculate the quantum corrections to 
polarization conversion and demonstrate that they can 
either increase of decrease conversion efficiency. 



A. Classical effects 

Under the condition Eq. ([55]) , the solution of Eq. ([T|) 
for f^{tf) reads 



n^) = % 






(36) 



The solution of Eq. ([2]) for s(ip) can be written in a 
straightforward way as well using Eqs. ([28|l . (|29p . In 
what follows we concentrate on the important case of ex- 
citation of the pure TE- or TM-state, i.e. 17^^ || g || x. 
Thus, 8° = gur/T), and 






9T_ 

V 



[2TrS{ip) + (l + nlT^)i^u] , (37) 



'■ flx^yT VU, 



Szif) = -l^^yTVU, 



where we introduced 

u = [1 + {VLrf - v{l 



riV 



m\ 



T 
Tl 



The angular distribution of the scattered polaritons is 
symmetric with respect to the rotation by an angle tt be- 
cause the longitudinal-transverse spitting contains even 
harmonics of the wave vector angle Lp^. Moreover, Sz 
pseudospin component appears to be proportional to the 
y-component of the splitting, Eq. ([57)1 . (see also Ref. 0), 
therefore the maxima of circular polarization appear at 
scattering angles equal to 7r/4, 37r/4, 57r/4 or Tit/ A. The 
sign of polarization is opposite in the adjacent maxima. 
Figure [3] shows the circular polarization degree of the po- 
laritons at the scattering by 7r/4. The dependence of the 
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FIG. 4: (Color online) The absolute value of the circular po- 
larization degree of polaritons as a function of Qt calculated 
for the different values of tq/ti. The scattering angle is tt/4. 
An inset shows a contour plot of the angular distribution of 
circular polarization degree, and arrow denotes the excitation 
point in the fe-space. 



polarization degree on fir is non- monotonous, its maxi- 
mum shifts to the lower values of f^r with decrease of the 
scattering time ti . The maximum value of the polariza- 
tion is observed for ti/tq — > oo, i.e. in the regime of a 
single scattering described in Rcf. [3. 

Multiple scattering leads to spin relaxation of polari- 
tons. At To ^ Ti, (ri^Ti)"^ one can see that the factor 
u reduces to 2/{Uti)'^ which is nothing but the ratio of 
the classical value of the transverse spin relaxation time 
T_L and the scattering time ti. The circular polarization 
degree in the scattered state 



P°(^) 



-—— = ny{ip)Ti— 
rm To 



is the smaller the shorter spin relaxation time. 



B. Quantum effects 

The quantum corrections to the particle and spin dis- 
tribution functions in the lowest order in Wq, W can be 
found similarly to Sec. IIVI The distribution functions 
are represented as the sum of zeroth order contributions 
and the first order corrections, see Eq. ([27|) . Functions 
f° and s° are given by Eqs. dM]), §7^. Thus, for p{ip) 
we have 

f\^)=TWo^[27T5{^-TT)-l] (38) 

The first term in the brackets describes the coherent 
backscattering, i.e. an increase by interference effects 
of the number of the particles scattered into the oppo- 
site from the source direction. The second term describes 
the coherent scattering by an arbitrary angle. The total 



The solution procedure for s^{ip) is analogous to the 
outlined above. We introduce the auxiliary function 



F\^) 






w 



s°{ip-n) 



(39) 



and the kinetic equation for s^ reduces to Eq. (|^5|) with 
F^ instead of F and s^ instead of s. Thus, s^ is given 
by Eq. (P5)l . One needs to average the solution over ip to 
obtain self-consistency equations for the components of 
s^ and then find s^((^). The result reads 



4(^) 



sliv) 



27rTW±S{ip 
4 "" ' 4~ 



w 



-Ao 



V 4 



y''' r< 



(40) 



1 gru^flyT 

z\i-; p 4 2 

Here the following quantities are introduced: 

A2 = AtW\\ {l-v)u- 2tW^p, 
Co = 8rVK|| (1 - v)v - AtWi_p, 



Bn 



SrWiiil - iy)iy + 8tW±{1 - v){v - I 



n\^) 



p = n^T^2 



4i^) + 2 + 4i^^ - 6iy. 



Quantum interference leads to the appearance of extra 
backscattered polaritons. The total number of backscat- 
tered particles is proportional to Wq, Eq. (|55|) . This in- 
crease is compensated by the coherent scattering by an 
arbitrary angle which leads to the decrease of the num- 
ber of scattered polaritons in all other directions than 
exactly backwards. An increase of the spin splitting at 
fixed Ti and tq leads to the decrease of Wo therefore the 
number of backscattered particles decreases. 

The coherent backscattering is also pronounced in the 
X pseudospin component. In this case the overall inten- 
sity of the backscattering peak is proportional to W± . 

The cross-linear polarization and the circular polariza- 
tion appear proportionally to ^x^y and Qy, therefore 
both of them vanish for ip = tt, i.e. for the detection in 
the backscattering geometry. However, the coherent scat- 
tering by an arbitrary angle leads to the modification of 
the polarization conversion efficiency. The circular polar- 
ization degree can be written as Pc{p>) = Pcif) + p]i{f)j 
where p" is the classical value of the circular polarization 
degree and the quantum correction p].{p) is given by 



number of particles is conserved, /„ f^{p)dp = 0. 



P^M=P°(^)(^-|J 



(41) 
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FIG. 5: (Color online) Quantum corrections to the circular 
polarization degree observed under optical spin Hall effect 
conditions. Solid curve shows the total correction, dash- 
dotted one shows first contribution, sl/s'^, in Eq. (|41|l . 
the dashed curve shows second term, f^/f", in Eq. (|41|l . 
To/n = 100. 



The relative value of the polarization conversion effi- 
ciency pI/Pc as a function of fir is plotted in Fig. El 

From Eq. (|4T|) it is clear that there are two contri- 
butions to the quantum correction to the polarization 
degree: first one arises due to the modification of spin 
distribution (s^) while second one is determined by the 
change of the number of particles in a given state. The 
latter correction is always positive because the coherent 
scattering by an arbitrary angle reduces the number of 
particles in a given state, see Eq. ([55)1 and dashed curve 
in Fig. [5] The former one can be either positive or nega- 
tive depending on the values of fir and tq/ti, see Fig. [S] 
For instance, if 17r <C 1 and tq » n both W\\ , W± are 
positive and Bq > 0- Therefore, quantum corrections in 
this regime increase the polarization as compared with 
the classical result. On the contrary, if fir ^ 1 both 
W\\ and W±{i' — 1 — O^r^) « —W±il'^T^ are negative 
and z-pseudospin component is decreased. Thus, inter- 
ference of polaritons can cither increase of decrease the 
polarization conversion efficiency. 

Physically, it can be interpreted as follows. The effi- 



ciency of the polarization conversion is strongly sensitive 
to the spin relaxation times. The correction to the trans- 
verse relaxation time can either be positive or negative 
depending on the sign of W\i and the value W±{i^T)'^, 
therefore the cc-pseudospin component can be preserved 
better or worse depending on the value of Qt. There- 
fore generated circular polarization and cross-linear po- 
larization may cither increase or decrease as a result of 
quantum interference. 

VI. CONCLUSIONS 

To summarize, we have studied in detail the spin- 
dependent quantum interference and classical multiple 
scattering effects in dynamics of exciton-polaritons. We 
have derived the quantum corrections to the collision 
integral of exciton-polaritons in the leading order in 
(fep/)"^. The quantum corrections are strongly sensitive 
to the value of the spin splitting of exciton-polaritons. 
Contrary to the case of electrons where the strong spin 
splitting can lead to the anti-localization, the quantum 
correction to the polariton diffusion coefficient is nega- 
tive. The quantum correction to spin z-component diffu- 
sion coefficient changes its sign from negative to positive 
with the increase of the spin splitting while the correc- 
tion to the diffusion coefficient of the in-plane spin com- 
ponents is negative. The relaxation of the longitudingal 
spin component is accelerated by the quantum interfer- 
ence effects and the relaxation rate of the transverse spin 
components can increase or decrease depending on the 
spin splitting value. The polarization conversion effi- 
ciency in the regime of the optical spin Hall effect can 
also be larger or smaller than the value predicted by the 
classical theory depending on the relations between the 
lifetime of polaritons, their scattering time and the value 
of the spin splitting. 
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